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ABSTRACT. For most applications in the control area, the 
standard practice is to approximate a nonlinear mathematical 
model by a linear system. Since the feedback linearizable 
systems contain linear systems as a subclass, we examine the 
procedure of approximating a nonlinear system by a feedback 
linearizable one. Because many physical plants (e.g. air- 
craft at the NASA Ames Research Center) have mathematical 
models which are "close" to feedback linearizable systems, 
such approximations are certainly justified. We introduce 
results and techniques for measuring the "gap" between the 
model and its "truncated linearizable part." The topic of 
pure feedback systems is important in our study. 

I. INTRODUCTION . In control design for nonlinear systems, 
the most common method is to approximate the nonlinear 
system by a linear system using the Taylor series trunca- 
tion. Thus, we approximate a nonlinear system by a linear 
one and design with respect to the linear system. Recent 
advances have shown that control design can be achieved 
using a much larger class (containing the linear ones) of 
systems, which are called "feedback linearizable." These 
are nonlinear systems which are feedback equivalent to 
controllable linear systems [1], [2], [3]. 
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Concerning applications, the feedback linearization 
approximation for totally automatic flight control is used 
in 14] and [5]. In [5] the particular case of the UH-lH 
helicopter is studied, and successful flight test results 
are discussed in [6] . 

The first author and R. Su [7] have examined the pro- 
cedure of approximating a general nonlinear system with a 
feedback linearizable one by introducing the concept of a 
pure feedback system. Every pure feedback system is feed- 
back linearizable, and for each feedback linearizable system 
there exists a state space coordinate system in which it is 
pure feedback (assuming generic controllability assump- 
tions) . These coordinates are called the s-coordinates and 
they are generated geometrically (see also [1] and [8]). 

In the s-coordinates we can quickly ascertain if the system 
is feedback equivalent to a linear system or not. If 
approximation is necessary, the "pure feedback part" is 
easily recognizable. 

Suppose we take a nonlinear system in its s-coordinates 
(1) s = f(s) + i f 1 u.g^s) 


with analytic vector fields on]R (or say an 

open set in R n containing the origin) , s = and u^,u 2 , 

...,u as controls. In general, the Taylor series approxi- 
mation of a single input system, m=l (we shall analyze 

. o 

multi-input also), has error 0( (Sj. ,s 2 , . . ■ ,s ) ) in terms 
of vector field differences. By O ( (s^ , s., , . . . , s ) 2 ) we mean 

no linear terms in these variables. However, the pure 

2 

feedback approximation has error 0( (s 3 ,s 4 , . . . , s ) ) . More- 
over, if g , [ f ,g] , . . . , (ad k f ,g) are involutive, 0<k<n-3, the 
pure feedback approximation error is O ( ( s^ , , . . . , s n _^) ) , 

and the vector fields in the original system and the pure 

feedback part agree when s„ = 0 . s, = 0, . . , , s .=0. If 

-• 4 n-K 

k = n-2, the system is pure feedback, and no approximation 
is necessary. 

In approximating a nonlinear system by its pure feed- 
back part, it is of interest to compare the state time 
responses of the system and its approximation. For this 
purpose we propose the Volterra series expansion of Fliess, 
Lamnabhi, and Lamnabhi-Lagarrigue [9]. 

An interesting class of partial differential equations 
consists of those failing to be elliptic, but possessing 
many desirable properties of elliptic equations. We con- 
sider the possibility of introducing the s-coordinates to 
find interesting coordinates for those hypoelliptic opera- 
tors studied in [10] and [11]. 

In keeping with the purpose of this Engineering Founda- 
tion Conference, we are presenting an overview of recent 
work and our thoughts and ideas about interesting problems 
and future directions. This concept is certainly reflected 
in the style and intent of this paper. 
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Sections II and III contain definitions, results, and 
examples concerning pure feedback approximation for the 
single input and multi-input systems, respectively. A 
discussion of geometrically generated coordinates for the 
study of partial differential equations is the topic of 
Section IV. Future research directions are mentioned in 
the final section. 

II. SINGLE INPUT SYSTEMS . We begin with a single input 
system 

(2) x = f (x) + ug(x) , 

where f and q are real analytic vector fields on some opeiu 
n . O 

set in Jr containing the origin. 

Dof ini f i at) 9 1 A flf thp foriTv 


f l (x l ,x 2 ) 
f 2 1 x 2 f x 3^ 


*n-l = f n-l (x l' x 2 X n> £ 

X n = f n (x l' x 2 X n> + u< ? n (x l' x 2 V * 

is called a pure feedback system . 

Let [f,g], (ad 2 f,g) = [ f , [f ,g] ] , . . . , (ad k f,g) = 

[f , (ad k_1 f ,g) ] denote Lie brackets involving the vector 
fields f and g. By L^y we mean the Lie derivative of a 

function y with respect to f; i.e. L^y = <dy,f>, <*,*> 

denoting the duality between one forms and gradients. We 
assume for the remainder of this section that g,[f,g],..., 

(ad n ^f,g) are linearly independent. 

For a pure feedback system (3) we define new state 
space coordinates ,y 2 , - . - ,Y n and a new’ control v by 


L f Y n-l 


v = L.v + uL V . 

n n 


We obtain the controllable linear system 


y, = y, 
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<S) 


Example 2.2 . On 


the system 
2 


( 6 ) 


4, = x, + x^ 

£ J 2 

. 2 
x 3 = x i + u 


Moreover, there is no 


is not a pure feedback system. 

3 

coordinate system on 3R in which this system appears as a 
pure feedback system. 

“ 3 is 


The usual Taylor series approximation about OeJC' 
x. = x 0 


(7) 


and the error between the vector fields in (6) and (7) is 


2 . . 2 , . 2 . 

, (x.) , (x, ) ) . 


Approximation by a pure feedback 


0( (x 3 
system 


( 8 ) 


yields an error difference in systems (6) and (8) of 
2 

0( (x 3 ) ). Moreover, the approximation of (6) by (8) is 

exact when x 3 = 0 . 


Example 2.3 . The system 


(9) 


X 2 + 


A' 


sin (x, -x_ ) 
l 3 


y n-l - y n 
y n = v. 

Hence system (3) is feedback linearizable or is feedback 
equivalent to the controllable linear system (5) . In 
general, for a feedback transformation we include non- 
singular state space coordinate changes (e.g. y^ ,y , . . . ,y n ) , 

additive state feedback (e.g. L f y n > , and nonsingular state 
dependent input space coordinate changes (e.g. uL^y n ) . 

In moving toward general results, the following two 
examples should be helpful. 
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2D 


s 

and, in 
Gi' 

s 

(2) 


with g. 

s 

find a < 

f 

i) 

J-Tk: 

-A 

ii) 


in [6] that near the origin there exists state space 
coordinates in which we do have a pure feedback system. 

Hence a pure feedback system is not invariant under 
coordinate changes on state space. Those nonlinear systems 
that can be reduced to pure feedback form can be classified. 


x = f (x) + ug (x) 

,g] , . . • . (ad n-1 f ,g) linearly independent, we shall 
ordinate system (called the s coordinates) so that 
if the system can be put in pure feedback form, it 
appears in this coordinate system, 
if the system cannot be put in pure feedback form, 
we approximate it by that part of the system in 
the s coordinates appearing in the form (3) . 

For every system (2) we have s coordinates and we can 
expand in a power series in these coordinates. 

Definition 2.4 . The pure feedback part of a nonlinear 
system is that part in the s coordinates which appears in 
the form (3) . 

We introduce the following coordinate system. Solve 
in order the following systems of ordinary differential 
equations with the indicated initial conditions. 


( 10 ) 


dx _ 
ds-^ 

(ad n 1 f,g) , x (0) ■ 

dx 

ds 2 

(ad n 2 f,g) , xls^ 

dx 

ds , 

= [f ,g] , x(s 1 ,s 2 , . 

n-1 

(X I' S 2 S n-2> 

dx 

ds 

<3 ' x ( s !' e 2 s n 

n 

x(s 1 ,s 2 , . . •»=„_!) • 

function 

theorem we invert 


x(s x ) 


n-2 ’ 


( 11 ) 


We set 


S 1 (X 1' X 2 V 

S 2 (X 1' X 2 V 


s n (x l' x 2 V • 


X 3 = (ad 11 L f,g) 


on B is not a pure feedback system. 


However, it is shown 
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X, = (ad ‘f,g) 


and define 


X n = g 


S k = {s = (Sj, s 2> . . . ,s n ) R : s m = 0, k + 1 s m s. n} 
for k = 1 ,2, . . . ,n. 

In the s coordinates, usual derivatives can be re- 
placed by derivatives with respect to the X^,j = 1,2,. ..,n. 

Theorem 2.5 [7] . In terms of the s coordinates the system 
(2) assumes the form 

n (s ) 1 

(14) § = f (0) + i | 1 (adV.f) s + ue n , 

where denotes restriction to S . . and e is an n 

I s j-1 3-1 V. 

vector whose only nonzero entry is a 1 in the n tn component. 

It is very easy to recognize the pure feedback part of 
system (14 ) . 

For a definition of an involutive set of vector fields 
and a statement of the Frobenius Theorem we refer the reader 
to [2]. 

Theorem 2.6 . If g, (f ,g] , . . . , (ad f,g) are involutive, k an 
integer, Osksn-3, then the vector field difference between 
system (14) and its pure feedback part is 
2 

0 ( (s 3 ,s 4 , . . . ,s k ) ). If k = n-2, then (14) is a pure 
feedback system. 

Proof. In the s coordinates, the manifolds S, are linear 

n k 

subspaces of R n . Moreover, on each the vector field 


takes for form X. 


, or in column vector notation 


Letting 


. th , 

. k place 


we find that (ad X n >f) = -rf,gj 
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= (ad f,g) implies f x 2 ' * ' " f n- (k+1) are 


independent of s 


n-k+1 ' " ’ ' 


Hence the error between system (14) and its pure feed- 
back part is clearly 0( (s 3 , s 4 , . . . , s n _ k > 2 ) . If k = n-2, we 

have a pure feedback system for (14) . !_] 

It is very important that a method be devised to com- 
pare the time responses of the state evolutions of system 

(14) and its pure feedback part. Given an initial condition 
s Q in the s coordinates and an input u we examine the 

difference between these time responses by applying a formal 
Volterra series expansion. 

We begin with system (2) and assume that we are in the 
s coordinates (i.e. x = s) . To 
(2) x = f (x) + ug (x) 

we add a real analytic output function y = h(x) and obtain 
the single input, single output system 
x = f(x) + ug(x) 

(15) y = h (x) . 

By the formula of Fliess et.al. [9] the formal 
Volterra series expansion of the system (15) is 


y(t) = J 0 L f hix c )^ + 


J. “ n ’ i * n 


^g^f n,x 0>- 


U(T 1 )dT 1 


I ' V ,. V - ' ' ~ 

I n l t li L L.i L, h(x„). 

->0 Jo^o ,v i ,v 2^^ f g f g f o 


(t-i,) (T,-r.) 


1 k /' 2 r 


— u ( i 2 )u(T 1 )di 2 di 1 


| 0 , .1 . , . k *0 L f°V‘f 1 • ’ • L g L f khiX 0 ) - 


a- !k )_ k ...(x i ) 0 

V-” V 


u(i k ) ulTj) dr k . . .dx^ 
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where x^eR is a point at which the system is defined. 

Our recommended method proceeds as follows: Take in 

order the outputs h = x^, h = x 2 ,..., h = x n , compute the 

Volterra series expansions for the system (15) and the 
system 

x = f(x) + ug(x) 

(17) y = h (x ) 

given by the pure feedback part, and compare the results 
for corresponding state time responses. 

We examine the effect of the involutivity assumptions 
of Theorem 2.6 on the Volterra series (again taking x = s) . 
From the proof of that result and the successive applica- 
tion of the formula 

L f <dh,g> = <dL f h,g> + <dh,[f,g]> 


we find that 


= L gL f h + <dh,[f,g]> 


L Lix. 
g f 1 

L L^x, 
g f 2 


0£v<k-l 


Lx,,, 0 

g k+1 

For any initial condition x„ and with h = x, , h = x_,... 

0 12 

= x k+3 in turn, the terms in the expansion (16) correspo 

ing to the Lie derivatives in (18) must vanish. 

Since the dynamical equation in (17) is pure feedback, 
it is feedback linearizable . This implies the set 

(g, If ,g] , . . . , (ad n ^f,g)r is involutive (see [12) and [2]) 
and 

L L^x, 0, 0<y<n-2 


Lx . 0 . 

g n-1 

As before, given a point x^ and h = x^, h = x 2> ...,h = 

in turn, the terms in the expansion (16) (with f replaced 
by f) corresponding to the Lie derivatives in (19) must 
vanish . 

The conditions in equation (19) common to those in (18) 
illustrate the importance of the involutivity assumptions 
in improving the errors in the corresponding Volterra 
series expansion differences for (15) and (17). Again, the 
outputs are taken successively to be h = x. , h = x 2 ,..., 

h = x . 
n 
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III. MULTI-INPUT SYSTEMS . We consider the nonlinear system 

(20) x = f (x) + ? u.g.(x) 

i=l 1 1 

with f,g.,...,g being real analytic vector fields on some 
open set in R n containing the origin, and g^ ,g 2 , . . . 

being linearly independent. Assume there exists a set of 
positive integers K x ' K 2 ' ' ‘ ' K m suc ^ that: 

i) the set C = (g^,(ad f.g.^), (ad f , g^ ) , 


(ad m f,g )} spans JR n on our open set containing 


ii) the span of C . 1C = the span of C.,j = 1,2, 

J k ,-2 ^ 2 

where = tg 1 ,...,(ad f .g^ ,y 2 , . . . , (ad J f,g 2 ) 

N4 '2 

(ad J f '5 m )} ' 

iii) ... (renumber g^ «g 2 » • • • ,9 m if 

necessary) . 

We fill a <jxm array by putting from top to bottom 

'■ i~i 

1) g^ ,( f ,g^ ),..., (ad f,gj) in the first column 

■ 2 -1 

2) g 2 , (f ,g 2 ) , . . . , (ad f,g 2 ) and 0's (if needed) 
in the second column 


< -1 

m) g m , [f ,g m ) , . . . , (ad f , g m ) and 0's (if needed) 

in the m*"* 1 column. 

.• -1 

Let X, - (ad f,g^j, the entry in the last row and first 
column, 

X 2 = the vector field entry in the last: row and 2nd 

column if it is nonzero, or the vector field entry 

th i — 1 

in the • ^-l 1 " row and 1st column if (ad f,g^) is 

the only nonzero entry in the last row. 


n in ^ 

Thus we start at the , row and first column and move 
x 

from left to right among the nonzero entries. Encountering 
a zero, we move up one row and return to the first column. 

The s coordinates are defined by solving m order the 
system of o.d.e.'s with initial conditions 
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dx 

di 1 X 1 

, X(0) = 0 


dx 

ds 2 x 2 

, x(s 1 ,0) = 

! X(s 1 ) 

^ = X 

ds n n 

, xfSjyS 

• • ,s n-l' 0) 

inverting 

(locally) to 

obtain 

s l (x l' x 2' 

...,x n ) 


s 2 (x 1 ,x 2 . 

' x n> 


S n < x i' X 2 ' 

...,x n ) . 



(22) j 

S n <X l’ X 2 V • 

The manifolds 

S Q = 0cR n 

= { s= ( s ^ , s 2 , . . . , s n ) JR n : s^ = 0, k+lsiin} 

are essential in the following result, which is a multi- 
input analogue to Theorem 2.5. 

Theorem 3 . 1 [ 7] . In terms of the s coordinates, the system 
(20) assumes the form 

(s . ) 1 

S = f(0) + -Ij j ^ 1 (ad 1 X j , f ) 

(24 ) ‘ 


+ k=n-m+l fe k + i£l jik+inT) (ad X j ,X k ) )u k-n+n 

j-1 

where e, is an n vector whose only nonzero entry is a 1 in 
the k component. 

We need a definition of pure feedback system for the 
multi- input case. Take the -:^xm array that we constructed 

above. For i.k=l ,2,.,., . , let n. = number of nonzero 

1 1 

elements In the (Kj-i+l)* row 

‘ k = n l +n 2 + ' " * +r! k ' 

For j = 2,3,...,<^ set 

y j = (x 1 ,x 2 , . . . ,Xg ) 

= (X 1 ,X 2 X n> 

Denote by g^ the j''* 1 component of the vector field g^. 

Definition 3.2 . The nonlinear system (20) is a pure feed - 
back system if it is of the form 

= f j <y 2 ) , 3 - 1,2, 

*j = f j (y 3^ ? * ? '1 +1 S 2 


i 
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f -i<Y„ >■ 


B _,+l, . . . ,n-ni 

<i i 


*j ' £ j<y< 1 ) + i £ l g ij (y i< 1 ' u i' i-n-a+1 

The block triangular systems found in [4] are a 
particular subset of the pure feedback systems. 

For a general nonlinear system (20) (perhaps not pure 
feedback) , we can move to the s coordinates as shown in 
equations (24) . 

Definition 3.3 . The pure feedback part of a nonlinear 
system (20) is that part xn the s coordinates which appears 
in the form (25) . 

Without proof we give a multi-input version of -Theorem 
2.6. The sets are defined in the assumptions following 

equation (20) . 

Theorem 3.4 . If each of the sets {g 1( (f .g^J , . . . , (ad x f ,g ) , 

9 2 <( f *9 2 ^ ' ‘ ' ‘ ' /9 2 ) ' g m' ' g m^' ' ’ ' ' ,g m^ an< * C j are 

involutive, where k is an integer, 0£k£K,-3, and j is any 
positive integer with Xj-2£k, then the vector field 

difference between system (24) and its pure feedback part 
2 

is 0( (Sg +1 ,...,Sp) ). Here p is the largest subscript on 


a nonzero vector field X i in the (k+1) 


row of our K^xm 


array. If k * <^-2, then (24) is a pure feedback system. 

We can compare the time state responses of the system 
(24) and its pure feedback part by using multi-input 
versions of the Volterra series expansions of [9). 

IV. PARTIAL DIFFERENTIAL EQUATIONS . Let L be a linear 
partial differential operator with real C°° coefficients and, 
for simplicity, assume L is second order. We also suppose 
that the principal part of L is the sum of squares of 
vector fields on an open set containing rhe origin inJR n . 

We let the principal part of L be f 2 + g? +...+g 2 , 
where . _ A m 

£ * a, (x) + a, (x) -r~- +. . .+ a (x) 

I ^x 2 ax. n dx 

1 i n 

g i = e ii (x, 3^r + s i2 (x) wi + --- + B in (x) dr 
• 12 n 


g m = 6 ml (x) 9x7 + e m 2 (x) 


t — + . . .+ 6 (x) 

*« 4«-4- n\4 dx» mn ox 

12 n 

and the squares mean that an operator is applied twice. 

In most studies the operator L is taken to be elliptic, 
but we are interested in the case that f can vanish on 
certain sets and ,<J 2 * - • • ,9^ are linearly independent with 

m<n. We assume the existence of integers k, k , sets 

12m 

C, and C. as in section III. 


C, and Cj as in section III. 
Definition 4.1 . An operator 


said to be hypoellipLic 


if Lu = f, where f is C on an open set U of IR n , implies 
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that u i« C“ on U. 

By the results of Hormander [10] (with extensions due 
to Rothschild and Stein [11] we have that our operator L 

with principal part f 2 +g 2 +...+g 2 is hypoelliptic since the 

vector fields in C are linearly independent. 

We remark that the s coordinates in Sections II and III 
for our real analytic systems of ordinary differential 
equations are also applicable for c“ systems. Thus we can 
view our operator L in the s coordinates as generated in 
section III. 

The relationship between controllability of systems of 
o.d.e.'s and hypoellipticity of p.d.e.'s has been well 
established in the literature. However, perhaps the 
special coordinates (e.g. the s coordinates) and the equi- 
valence results from o.d.e.'s have not been applied to 
yield nice coordinates and equivalence criteria for opera- 
tors like L. For example, the theory of Krener [8] and 
Respondek [13] for state space equivalence of systems when 
applied to p.d.e.'s produces the following theorem in the 
m=l (with g=g^) case. 

Theorem 4.2 . For the second order partial differential 
operator with principal part f 2 +g 2 on K n , there exist a 

(local) coordinate system on in which the principal part 
2 2 

appears as (Ax) +b , where 

A is a real constant n by n matrix 
b is a constant vector 

if and only if [(ad r f,g), (ad s f,g)] = 0 for (Kr,s<n. Here 

g ;; 

Ax and b denote the coefficients of -r— , . 

X 1 x 2 n 

In fact, under the assumptions of the theorem, A can 
be in rational canonical form and 



Perhaps we can find "geometrically interesting canonical 
forms" for the class of partial differential equations 
under consideration. 

V. FUTURE DIRECTIONS . The s coordinates of section II and 
III are generated by solving systems of ordinary differen- 
tial equations. For feedback linearizable systems, 
symbolic and numerical methods exist in some cases for con- 
structing the feedback linearizing transformations [4], 

[14], the major task being to find the s coordinates. G. 
Blankenship et. al. at the University of Maryland are 
developing an expert system to construct such transforma- 
tions. The process of finding the s coordinates for general 
systems seems quite difficult. However, for many physical 
systems the mathematical model is often "near" pure 
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feedback form and "close" to being in the s coordinates. 

In this case, a "good" initial guess for the s coordinates 
is easily obtained and an appropriate iteration scheme is 
recommended. Efforts to handle more difficult nonlinear 
systems are presently underway by G. Meyer of the NASA Ames 
Research Center. 

When an approximation method is used for nonlinear 
systems of ordinary differential equations, the ultimate 
goal is to achieve close time responses between the 
original system and the approximating system. The Volterra 
series expansions mentioned in Section III seem ideal. We 
shall develop a symbolic manipulation program to generate 
and compare Volterra series expansions for nonlinear con- 
trol systems and their pure feedback approximations. 

An important numerical method in the engineering and 
mathematical approach to controlled systems of partial 
differential equations is the finite element method. How 
do the geometric hypoellipicity conditions of Section IV 
influence the use of finite elements? Moreover, one often 
does have a system of partial differential equations to 
model a physical system, but a computer generated finite 
element model. Is it possible to find the desired 
geometric conditions in these finite element models? 
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